Conductance of a hydrogen molecule 
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A recent experiment reported that a single hydrogen molecule can form a bridge between Pt 
electrodes, which has a conductance close to one quantum unit, carried by a single channel. We 
present density functional calculations explaining these experimental findings. We show that the 
symmetry of the molecular orbitals selects a single conduction channel. The transmission of this 
channel is close to unity due to a combination of charge transfer between H2 and the Pt contacts 
and the strong hybridization between the bonding state of the molecule and the d band of the Pt 
leads. 
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The recent advances in nanofabrication have triggered 
the hope that electronic devices can be shrunk down 
to the single molecule scale ||. Single molecules 
have been shown to perform functions analogous to 
those of the key microelectronics components such as 
switches ||, |[ ||, rectifiers Q and electronic mixers 0. 
However, the future of molecular electronics depends cru- 
cially on our understanding of the transport mechanisms 
in single- molecule junctions. From the experimental side 
an unambiguous contact to a single molecule is difficult 
to achieve ||, and in many cases the measurements are 
not reproducible. On the other hand, the discrepancies 
between theory and experiment and even between differ- 
ent theories are notorious 

From this perspective the measurement of the conduc- 
tance of an individual hydrogen molecule reported by 
Smit et al provides a valuable opportunity to an- 
alyze the emerging concepts on the electrical conduction 
in single-molecule devices in the perhaps simplest pos- 
sible system. In Ref. [[lo| it was shown that a single 
hydrogen molecule can form a stable molecular bridge 
between platinum contacts. In contrast to results for or- 
ganic molecules [|nj |l3| , this bridge has a conductance 
close to one quantum unit, carried by a single channel. 
This result belies the conventional wisdom because the 
closed-shell configuration of H2 results in a huge gap be- 
tween its bonding and antibonding states, making it a 
perfect candidate for an insulating molecule. Thus, the 
observations of Ref. jl(J pose an interesting theoretical 
challenge and their understanding can help to elucidate 
some of the basic transport mechanisms at the single 
molecule scale. 

In this paper, we develop (i) a simple theoretical model 
to understand the conduction mechanism through the hy- 
drogen bridge which (ii) we match to quantitative density 
functional theory (DFT) calculations of the conductance. 
The interplay of these two approaches permits us to iden- 
tify the two relevant mechanisms that give rise to the 
large conductance of the hydrogen molecule. First, the 
catalytic activity of platinum is responsible for a signif- 
icant charge transfer between H2 and the Pt contacts, 
which moves the bonding state of the molecule towards 



the Fermi energy. Then, the strong hybridization with 
the d band of the Pt leads provides a large broadening 
of the bonding state, finally allowing for a high trans- 
mission. These two key findings are at variance with the 
central conclusion of the theoretical analysis presented in 
the experimental paper. Additionally, we show that, due 
to the symmetry of its molecular orbitals, H2 filters out 
only one of the Pt conduction channels. We discuss the 
fundamental implications of these findings for the proper 
description of the transport properties of single-molecule 
devices. 

In order to investigate the transport through the 
hydrogen bridge reported in Ref. [[[(J we have per- 
formed DFT calculations of the linear conductance us- 
ing the method described in Ref. [Q, similar in spirit 
to Refs. |H| In this approach, we first decom- 

pose the Hamiltonian of a molecular junction as H = 
Hl + Hr + Hc + V , where He describes the "central clus- 
ter" of the system, which contains the molecule plus part 
of the leads, Hl,r describe the left and right electrode 
respectively, and V gives the coupling between the elec- 
trodes and the central cluster. The electronic structure of 
the central cluster is calculated within the DFT approxi- 
mation |l7j . The left and right reservoirs are modeled as 
two perfect bulk crystals of the corresponding metal us- 
ing the tight-binding parameterization of Ref. Q . The 
inclusion of part of the leads in the ab initio calcula- 
tion assures the correct description of (i) the molecule- 
electrodes coupling, (ii) the charge transfer between the 
molecule and the leads and (iii) the lineup of the molecu- 
lar levels relative to the metal Fermi level [|l9| . The Fermi 
level is set naturally by the highest occupied molecular 
orbital (HOMO) for a sufficiently large number of metal- 
lic atoms in the central cluster. 

As in the experiment we concentrate on the anal- 
ysis of the linear conductance, which is given by the 
Landauer formula G = GoTr{tft} = Go^2 i Ti, where 
Go = 2e 2 /h is the quantum of conductance, i is the trans- 
mission matrix and the TVs are the transmission eigenval- 
ues at the Fermi energy Ep. The transmission matrix is 
given by i(e) = 2 f 1 ^ 2 (e)G r c (e)t]( 2 (e) with the scattering 
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rate matrices Tl,r = ^- m (^L,R)- Here Yil,r are the self- 
energies which contain the information of the electronic 
structure of the leads and their coupling to the central 

cluster. Finally, the Green functions of the central clus- 
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ter are given by (?c(e) = el — He — ^z,(e) — ^(e) 

For the case of non-orthogonal basis sets it is straightfor- 
ward to implement the overlap matrix into these formulas 
via a Lowdin transformation (see e.g. Ref. ]20|). 

In Ref. (l0| a break junction was used to produce pure 
Pt contacts of atomic size. Then, it was shown that the 
transport through these contacts was markedly altered in 
the presence of hydrogen. Following the experiment we 
first analyze the transport through a pure Pt single-atom 
contact. In the DFT calculations we use the BP86 func- 
tional and the basis set of Christiansen et al. [p2| . 
For the description of the Pt reservoirs we use a basis 
with the 5<i, 6s, 6p orbitals. Fig. p] shows the transmis- 
sion and local density of states (LDOS) Q projected 
into the central atom of the structure shown in the in- 
set. We find that the HOMO of the central cluster is at 
-5.44 eV, which is close to the Pt work function ~ -5.64 
eV. The gap between the HOMO and the lowest unoc- 
cupied molecular orbital is almost closed (0.12 eV). The 
total transmission at the Fermi energy is T to t = 1-47 in 
agreement with the conductance histogram of Ref. [ [l0[ , 
which shows a large peak at 1.4-1.8Go. This transmis- 
sion is made up of five conduction channels, which is due 
to the contribution of the Pt d orbitals. This is again a 
manifestation of the fact that in single-atom contacts the 
number of channels is controlled by the valence orbitals 
of the central atom [^4], p5| . In particular, the dominant 
channel is a combination of the s and d z i orbitals p6[ . 

Let us now study how the presence of H2 modifies the 
conduction. Usually the lack of knowledge of the pre- 
cise geometry of the molecular junction complicates the 
comparison between theory and experiment. However, in 
Ref. fuj the presence of H2 was identified by means of 
the signature of its vibrational modes in the conductance. 
This information establishes clear constraints on the ge- 
ometries realized in the experiment. In this sense, we 
only consider configurations which are compatible with 
the observed vibrational modes. The most probable con- 
figuration is shown in the inset of Fig. ||, where the H2 is 
coupled to a single Pt atom on either side (top position). 
In this geometry the vibrational mode of the center of 
mass motion of H2, which is the one seen in the exper- 
iment, has an energy of 55.6 meV, lying in the range 
of the experimental values. In Fig. |] we also show the 
transmission and the LDOS projected into the orbitals of 
one of the H atoms. The total transmission at the Fermi 
energy is T tot = 0.86 and it is largely dominated by a sin- 
gle channel, in agreement with experimental results. We 
would like to draw the attention to the following two fea- 
tures in the LDOS: (i) the bonding state of the molecule 
appears as a peak at ~ 6 eV below Ep and the anti- 
bonding state, not shown in Fig. ||, is located at ~ 18 eV 
above Ep. This indicates that the molecular character 




FIG. 1: Transmission and LDOS at the central atom as a 
function of energy for the pure Pt structure shown in the 
inset. The LDOS is projected into the Pt valence orbitals. 
The Fermi energy is set to zero. Inset: Pt one- atom contact 
consisting of a fee structure with bulk interatomic distances 
grown along the (111) direction. The blue atoms represent 
the bulk atoms used to model the reservoirs. 



of H2 is largely conserved, (ii) Around the Fermi energy 
the gap between the molecular states is filled due to the 
strong hybridization with the Pt leads. 

We now proceed to explain the underlying physics with 
the help of a simple model. To illustrate the ingredients 
necessary to reproduce the full DFT results we shall pro- 
gressively sophisticate the model. We describe H2 with 
a two-sites tight-binding model, see Fig. |(a). In this 
scheme eo represents the Is energy level of H and tn is 
the hopping connecting the H atoms. This hopping is 
simply related to the splitting between the bonding (e+) 
and the antibonding state (e_) of the molecule, namely 
e± = eo ± tn, and its value is ~ —12 eV. The molecule is 
coupled symmetrically to the leads with a single hopping 
t. Obviously, within this model the conductance is made 
up of a single channel, but let us postpone the discussion 
of this point for the moment. The transmission is given 

by 



T(e) 



4TH 2 H 



[( e _g + )2 +r 2] [{e _~ e _ ) 2 + T 2y 



(1) 



Here, e± = eo ± tn — t 2 Ke{g a } are the renormalized 
molecular levels, g a (e) being the advanced Green function 
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FIG. 2: Transmission and LDOS projected into one of the H 
atoms as a function of energy for the Pt-H2-Pt structure, the 
central cluster of which is shown in the inset. At Ef T\ — 0.83 
and T2 = 0.03. The H-H and Pt-H distances are 0.8 A and 
2.1 A respectively. We use the cc-pVDZ basis set for H. 
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FIG. 3: (a) Schematic representation of the model described 
in the text, (b-c) Transmission and DOS projected into the 
bonding (solid lines) and antibonding (dashed lines) states 
vs. energy for different values of the scattering rate T. (d) 
Bulk DOS of the Pt reservoirs, (e-f) Transmission and LDOS 
projected into one of the H atoms for two values of the hop- 
ping £, assuming that the leads are modeled by Pt bulk atoms, 
tu =-12 eV and eo = 6 eV. In panels (d-f) the Fermi energy 
is set to zero. 



which describes the local electronic structure of the leads. 
The scattering rate T, which determines the broadening 
of the molecular levels, is given by T(e) = t 2 Im{g a } = 
7rt 2 p, where p(e) is the LDOS of the Pt contacts. Let 
us first assume that T is energy independent and that 
the levels are not renormalized (e± = e±). In Fig. |](b) 
we show the transmission as a function of energy for 
different values of T in units of tu- We also show in 
Fig. H(c) the corresponding DOS projected into the bond- 
ing and antibonding states of H2, which are given by 
p± = r/7r{(e — e±) 2 + T 2 } respectively. Taking into ac- 
count the huge value of tn , one might naively expect the 
curve for T = O.Obtn to represent the relevant situation. 
Assuming that H2 remains neutral, the bonding state is 
occupied by two electrons and Ef = eo- In this sim- 
ple picture H2 would be insulating, in clear disagreement 
with the observations of Ref. |10[| . Our DFT calcula- 
tions indicate that there are two main ingredients miss- 
ing in this simple argument. First, H forms a covalent 
bond with Pt by sharing electrons. The DFT calculations 
show that every H atom donates ~ 0.1 2e to Pt. This im- 
plies that the Fermi level lies closer to the bonding state. 
With this charge transfer the transmission raises signifi- 
cantly, but it is not yet sufficient to reproduce the DFT 
results. Thus, a large broadening (T) of the bonding 
state is still needed. As suggested by the DFT calcula- 
tions, this is provided by the good Pt-H2 coupling and 



the large DOS around the Fermi energy coming from the 
d band of Pt. We test this idea assuming that g a is the 
bulk Green function of Pt. The Pt bulk DOS is shown 
in Fig. ||(d). We also show in Fig. ||(e-f) the transmis- 
sion and the LDOS projected into one of the H atoms 
for two values of the coupling to the leads t. One can see 
that for realistic values of t « 1-2 eV, the transmission at 
Ef can now indeed reach values close to one. Therefore, 
we conclude that the high conductance of H2 is due to 
the charge transfer between H2 and the Pt leads and the 
strong hybridization between the bonding state and the d 
band of Pt. This mechanism is not exclusive of Pt and 
it must also operate in other transition metals, as it was 
shown experimentally for Pd (see Ref. flC[]). 

Let us now address why only a single channel is ob- 
served. In view of the model described above a simple ex- 
planation could be that all the current flows through the 
Is orbital of the closest H atom to the Pt atom. However, 
the DFT calculations show that while the Pt-Pt coupling 
is negligible, this is not the case for the Pt and the second 
H atom. In principle there are two paths for the current, 
i.e. two channels may occur. Due to the spherical sym- 
metry of the H orbitals there is only coupling to the s 
and d z 2 orbitals of Pt. In addition, the Hamiltonian ma- 
trix elements between these orbitals fulfill approximately 
the condition Hi s H2,d 2 = 2 ^2,sj where 1 and 2 
stand for the first and second atom of H2. Therefore, the 
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FIG. 4: Transmission as a function of energy for the structure 
shown in the inset. The H-H and Pt-H distances are 0.8 A 
and 1.86 A respectively. The energy of the vibrational mode 
of the center of mass motion of H2 is 65 meV. 



coupling matrix between the H2 and the Pt leads is sin- 
gular, which implies that one of the transmission eigen- 
values vanishes. This condition of the hopping matrix 
elements effectively reduces the rank of the transmission 
matrix. In other words, one of the molecular states does 
not couple to the metallic states due to symmetry reasons, 
reducing the actual number of channels to one. This phe- 
nomenon of symmetry-inhibited modes is well-known in 
the contexts of light propagation in photonic crystals 
and sound propagation in periodic structures j p8[ . This 
symmetry induced destructive interference is a unique, 
but generic feature of the top position and thus a direct 



indicator of the relative orientation of the H2 axis to the 
electrodes. 

In principle there are other geometries compatible with 
the vibrational modes analysis. However, based on the 
channel analysis performed in the experiment many ge- 
ometries can be ruled out. As an example we consider 
the situation sketched in the inset of Fig. |4| where each H 
is bound to three Pt atoms (hollow position). Indeed this 
configuration is suggested in studies of the chemisorption 
of H on Pt surfaces [p9[ . The conductance is carried by 
up to 7 individual channels (see Fig. Q). Due to the short 
distance between the Pt leads most of the current is by- 
passing the H2 going directly from Pt to Pt. This analysis 
allows us to conclude that this type of geometries is not 
realized in the experiment and illustrates the importance 
of the channel analysis. 

In conclusion, we have presented a theory for the con- 
ductance of a hydrogen bridge between Pt contacts ex- 
plaining the experimental observations of Ref. ]Iq[ . We 
have shown how H2 filters out one of the conduction chan- 
nels of Pt with nearly perfect transmission. Our anal- 
ysis of this ideal test system illustrates that ingredients 
such as charge transfer and the electronic structure of the 
metallic contacts are essential for the proper description 
of the electrical conduction in single-molecule devices. 
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